RARE DECAYS OF HEAVY QUARKS SEARCHING 
GROUND FOR NEW PHYSICSQ 



PAUL SINGER 

Department of Physics, Technion-Israel Institute of Technology, Haifa, 

Israel 

The search for new physics beyond the standard model is proceeding 
nowadays intensively along experimental and theoretical lines. We review 
here the sector of charm radiative decays in this context. The calculation 
oi D ^ V'j, D transitions reveals their unequivocal dominance 

by long-distance contributions. On the other hand, the beauty-conserving 
charm-changing electroweak transition 5*7 is shown to have unique 

properties which make it a promising avenue in the search for new physics. 
We describe a calculation of short- and long-distance contributions to this 
decay which finds them to be of comparable size. The branching ratio of this 
decay in the standard model is estimated to be ~ lO"*^. 



PACS Numbers: 12.39.Fe; 12.39.Hg; 13.25.-k; 13.30 Hq. 
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1. Introduction 

The standard model (SM) of strong and electroweak interactions [1], 
based on local gauge invariance with respect to the gauge group SU(3)c x 
SU(2)^ X U(l)y, is presently in excellent shape. All experimental data are in 
agreement with SM and its relentless success is providing an ever increasing 
challenge to both experimentalists and theorists. The only missing fixture 
of the model is the SM Higgs boson, for which the existing experimental 
searches put a lower bound on its mass of 88.6 GeV/c^ [2]. 

Despite its remarkable success, it is generally believed that the SM is 
in fact an effective theory at the energies presently accessible. This belief 
is fueled by the fact that the SM has about 20 arbitrary parameters and, 
moreover, there is no satisfactory explanation for many of its salient features. 
For instance, why does the gauge group of interactions have the structure 
expressed by its three factors? Why three generations of fermions? Why 
left-right asymmetry? How to explain the observed spectrum of quark and 
lepton masses and the pattern of mixing angles? 

This situation has led theorists to propose many paths for the possible 
extension of the standard model. I do not plan to go into any detail here on 
the variety of possibilities, which was reviewed at many conferences [3,4] , and 
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I shall restrict myself to the mention of a few of the more widely-discussed 
proposals: supersymmetry [5], especially the "low-energy" Minimal Super- 
symmetric Standard Model [6] which is considered as a most likely possibility 
at the Fermi scale, grand- unified theories [7], right- left models ]8], two Higgs 
doublet models [9], flavour changing neutral Higgs models [10], multiple Z° 
bosons [11], and anomalous triple gauge boson couplings [12]. 

An important feature of the standard model is the flavour symmetry, 
as the gauge interactions do not distinguish among the three generations of 
leptons and quarks. In practice, this symmetry is broken, as it is evident from 
the pattern of masses and mixing angles of the SM fermions. The Higgs boson 
is an agent of flavour symmetry breaking in SM via its Yukawa couplings to 
fermions. However, this flavour symmetry breaking is realized in a particular 
way; the tree-level neutral couphngs of the Higgs boson, as well as those of 
the photon and of the Z° boson are all flavour diagonal. The observed neutral 
flavour-changing processes on the other hand are rather small, being made 
possible in SM by loop graphs only; as such their magnitude is determined by 
the values of quark masses and of the CKM matrix. In view of the smallness 
of flavour-changing neutral-current (FCNC) in the standard model, FCNC 
transitions are usually considered to be a fertile ground for the search of 
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processes induced by new physics, which does not automatically suppress 
such processes. The charm sector plays a special role in this respect, since 
as a result of the effectiveness of the GIM mechanism in this sector, the 
short distance SM contributions to certain charm processes are very small. 
Accordingly, D° — D° mixing and rare charm decays have been singled out 
[4,13] as attractive candidates for the discovery of new physics effects. 

In this lecture we consider the potential of the electroweak penguin tran- 
sitions c — > in the search for new physics. To begin I shall review shortly 
the SM physics of the single quark transition Q — > 57, then I shall present the 
status of short- distance and long-distance contributions in processes driven 
by c — > W7 and c — > ut^l" transitions and finally I shall describe our recent 
work [14] which singles out the —>■ 5*7 decay as a unique tool for the 
search of effects beyond the standard model. 



2. The flavour changing Q ^ transition 

Flavour-changing photon transitions from a heavy quark Q to a light 
quark q are induced by loop diagrams and are a basic feature of the standard 
model [15], generally recognized as "electroweak penguins". Typical transi- 
tions are s — > dj, b —>■ s{d)j for which up-quarks contribute in the loop, and 
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c — > uj, t — > c(m)7 which are driven by down-quarks in the loop. 

The amphtude for such transitions, with the quarks Q^qon. the mass-shell 
is given by [15] 



+ F2j{k )ia^yk Mq— + TTlq 



u{Q) . (1) 



Fi , F2 are the charge-radius and magnetic form factors respectively and 
Vab are CKM matrices; Fi,F2 were first calculated in the electroweak SM 
by Inami and Lim [16]. The Fi term does not contribute to decays with 
real photons, however, it is relevant in leptonic decays hke B — > X{s)£i, 
K —>■ ttH, D — > VU. In order to compare the calculations of these processes 
with experiment, one must complement the electroweak SM calculation by 
the inclusion of QCD corrections [17]. In this section we shall mention the 
s — > (^7 and 6 ^ 57 transitions and in the next section we turn to the charm 
sector in more detail. 

The contribution of s ^ ^7 to various radiative X-decays [18,19] and 
hyperon decays [20-25] has been studied extensively in the last twenty years. 
As it turns out [26] radiative processes which are in the ~ (10^'^ — 10~^) 
range of branching ratios like — > 7r"'"7r°7 , — > 7r"'"e'*"e~, — > p7, 
5~— >E~7 have both short-distance and long-distance contributions and the 



latter are dominant; this prevents a direct and trustworthy check of the SM 
or of deviations from it in these decays. In order to investigate the short 
distance s d"j" transition one must turn to very rare decays [27], hke 
K'^ n'^uP, K1 7r°e'^e~, — > n^uP. In these, the short-distance 
contribution is prominent and the QCD corrections to the decay amphtudes 
have been estimated [28]. The most frequent of these is K'^ — > n'^uP, which 
is expected [28] in SM with a branching ratio Br(ir+ -k+vP) = (9.1 ±3.8) x 
10^^^. Recently [29], one event has been detected in this channel, which gives 
Br(ir+ 7r+i/p)exp = (4.213:5) x 10~^°. The other two decays are expected 
with branching ratios of the order of 10~^^ and one must wait for the planned 
experiments in order to find out whether the s — > cT' 7" and the box diagrams 
involved of SM give an accurate picture for these transitions. In the domain 
of hyperon radiative decays a similar situation prevails [25] ; however, there 
might be an exception as it appears [22] that the yet unobserved H~7 
decay is affected in a measurable manner [24,26] by the SM single quark 
s ^ d'j transition. 

Although the s ^ dj was the first to be investigated with the aim of 
relating it to the observed radiative decays of kaons and hyperons, it is the 
6 — > 57 transition [30] which has been the center of attention during the last 
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dozen years. Since it was pointed out [31] that the enhancement provided 
by QCD corrections to 6 — > 57 (in which the top quark in the loop gives 
the main contribution) would bring the inclusive B — > Xgj and exclusive 
B ^ K*^ decays into the realm of observability, a considerable amount of 
theoretical activity has proceeded alongside the experimental observation. 
The CLEO collaboration was the first to measure the inclusive rate [32] 
Br{B Xs-f) = (2.32 ±0.57± 0.35) x 10"^ as well the exclusive (charged and 
neutral) decay [33] Br(5 K*'y) = (4.5 ± 1.5 ± 0.9) x 10"^ The theoretical 
effort has been directed on the one hand towards a best determination of the 
QCD corrections to the inclusive process in SM and on the other hand to 
establishing the limitations imposed by the observed rate on various "beyond 
the standard model" theories. For typical Refs. on the latter effort see [34]. 
The latest theoretical calculations within the SM[35] give Br(S — > ^57) = 
(3.32 ±0.30) X 10~^ which should be compared with two recent experimental 
results: the CLEO update giving [36] Br (6 57) = [3.15 ± 0.35 (stat) ± 
0.32 (syst) ±0.26 (mod)] x 10~^ which is derived from an analysis of 3.3 x 10^ 
BB pairs andthe ALEPH result [37] of [3.11±0.80 (stat) ±0.72 (stat)]xlO-l 
Obviously, the agreement with the SM is impressive. 

There are two remarks to be made here. Firstly, the conclusion on the ex- 
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cellent agreement with SM assumes that LD contributions are small, which is 
indeed the result of many calculations (approximately 5-10%) [24]. Secondly, 
we await for experimental results on the complementary process b si~i~^ 
(including B K*£+e-, B Kl+l') which should be compared with SM 
theoretical expectations of a branching ratio in the 10~^ range. 

Before turning to the charm sector, we conclude that the study of the 
Q ~^ 0.1-1 Q ~^ qt^i' transitions in SM is waiting for the measurement of 
very rare decays in the strangeness domain, while in the beauty sector, where 
experiments are available, the standard model does very well so far. 

3. Short distance c — > ti7 and c — > ut^l~ 

The c ^ w) transition is induced by the electroweak penguin with the 
down quarks running in the loop. In the absence of QCD corrections this 
transition is extremely small as a result of the small masses of the quarks 
in the loop and the smallness of the CKM factors. The electroweak SM 
calculation [38] gives for this strongly GIM suppressed transition a branching 
ratio of ~ 10^^^ only. Including the QCD corrections at the leading log 
approximation [38], the C7 Wilson coefficient of the a^j, operator gets the 
admixture of Ci, C2 Wilson coefficients and the amplitude is increased by two 
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orders of magnitude, giving a branching ratio of about 10~^^. The calculation 
of the complete two-loop QCD corrections [39] leads, after using unitarity of 
CKM, to the folowing effective Lagrangian 

C7(mc) = 0.0068 - 0.020i . (2) 

^Prom this expression, another increase of two orders of magnitude in the 
SD amplitude obtains, giving rise to ^^r{c u^)/V{D°) ~ 2.5 x 10"^ This 
implies that in exclusive modes, like D Vj, the SD contribution to the 
branching ratio would be about (3 — 5) x 10~^. In order to ascertain the 
possibility of detecting the SD transition, one must now consider the size of 
the LD contribution. 

The SD amplitude for c — > ut^l~ can be obtained from the general 
electroweak amplitude [16], and the explicit expression for the effective La- 
grangian after certain simplifications is [40] 

A = -0.065 . (3) 

This electroweak transition is not strongly suppressed, in contrast to c — > 1^7 
and although the QCD corrections have not been evaluated explicitly, they 
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are not expected to change the value of A appreciably [40]. From (3) one 
finds s°r(c ue+e-)/r{D'') ~ 3 x lO-^. Hence, hke in the c ^ u-y case, 
one has to ascertain the LD contribution before one may use these leptonic 
decays for checking the standard model. 



4. D-mesons radiative decays - the long distance aspect 

Several treatments have addressed recently the problem of estimating 
LD contributions to radiative D decays. These approaches include a pole 
model [38], a quark model [41] , the use of QCD sum rules [42] and effective 
Lagrangians [43]. Already from these works one learns that the D 
decays are expected to have branching ratios of the order of 10^"^ — 10^^, 
much larger than from the SD part. 

In a more comprehensive and systematic treatment for these decays [44] 
we used an effective hybrid Lagrangian combining heavy quark symmetries 
and chiral symmetry [45] to calculate nine decay modes of the D Vj type. 
The effective nonleptonic Lagrangian used is given by 

G 

Lld = -^Vuqy*q^[ai{uqiY{qjc)^ + a2{uc)^{qjqiY] (4) 

and for the QCD-induced constants 01,02 we take Oi = 1.26, 02 = —0.55 as 
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determined [46] from nonleptonic D decays. In order to evaluate the matrix 
elements of (4) we use the factorization approximation for the (yVo\{qiqjY{qkc) ij\D) 
amplitudes. 

The general gauge invariant amplitude for the decay D{p) — > V{pv)+'^{k) 



IS 



+ i [i^lv) ■ • P(y)) - {P(V} ■ k) } Apy . (5) 

In Ref. [44] all diagrams contributing to Apc, Apv are classified and their 
explicit expressions are presented. In Table 1 below we give the predicted 
widths [44] as well as the existing experimental upper hmits [47]. Since 
the amplitudes contain several terms, with unknown relative phases, we can 
present only their expected range. The first two decays in the Table are 
Cabibbo-allowed, the next five are Cabibbo-forbidden and the last two are 
doubly forbidden. To give an indication, the photon energy in the first two 
decays is 717 and 834 MeV respectively. As it is obvious from Table 1, in all 
these decays the LD contribution masks totally the SD one - preventing the 
detection of deviations from it by orders of magnitude. 
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D ^V'-^ Transition Br Ratio x 10^ [44] Exp. limits [47] 



D" K*° 6-36 < 7.6 x 10"^ 

Df p+ 20-80 

p° on < 2.4 X 10-^ 

D^^u) 0.1-0.9 < 2.4x10-^ 

D° ^ip 0.4-1.9 < 1.9 X 10-^ 

D+ p+ 0.4-6.3 

Dt K*+ 1.2-5.1 

D+ K*+ 0.03-0.44 
D° K*" 0.03-0.2 



Table 1 

Turning now to decays of type D — > Vi~^i~, these were also calculated 
recently [40] using generally the same theoretical framework [45] as for D ^ 
Vj transitions. Since the SD transition (Eq. 3) is considerably larger here 
than in the c — > it7 case before the application of the QCD corrections, one 
could expect that the gap between SD and LD contributions is narrower for 
the leptonic decays in the SM. Such a situation could open the window to 
new physics. 

The authors of Ref. [40] have considered the same hadronic transitions 
as in Table 1. The SD contribution due to c — > ul^l~ is present in the five 
Cabibbo suppressed decays — »■ a;°, v?'')^^^^, p+7, K*^^ 

while in the other four decays signals for new physics might come from more 
exotic contributions. The calculation is performed [40] again using factoriza- 
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tion for matrix elements of (4) which leads to three classes of diagrams: the 
annihilation contribution, the V^-spectator part and the ^-spectator part, 
where V is the final state particle and an intermediate vector meson 
{p,uj,ip). There are thus two kinds of LD contributions: the resonant mech- 
anism, where in addition to V also Vo is produced in the final state and 
converts to a photon through vector meson dominance, and a nonresonant 
mechanism with the photon emitted directly from the initial D state, as pre- 
scribed by the structure of the hybrid lagrangian [45]. The latter should 
contain in our approach also possible contributions from intermediate cc 
states. The predicted branching ratios for D — > V , including SD -|- 
LD contributions, and the exisiting experimental upper limits are given in 
Table 2. The range in column two is due to coupling parameter uncertainties. 



D Vn+n' Calculation [4] of Br(LD+SD) Exp, hmits [48] 



1.6 - 1.9) X 10-6 
3.0 - 3.3) X IQ-^ 



< 1.18 X 10-3 



Df K*+ 



3.5 - 4.7) X 10-^ 

3.3-4.5) X 10-^ 
6.5-9.0) X 10-^ 



< 2.3 X 10-4 

< 8.3 X 10-^ 

< 4.1 X 10-4 

< 5.6 X 10-4 

< 1.4 X 10-3 




5.0 - 7.0) X 10-^ 



D+ K*+ 



3.1 - 3.7) X 10-^ 
4.4 - 5.1) X 10-9 



< 8.5 X 10-4 



Table 2 
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The short-distance contributions alone are ~ 10 ^ for D° —>■ p°{uj°)ii'^ fi , 
5 X 10-9 for D+ P'^IJ.'^I^' and 1.6 x 10"^ for D+ K*+fj,+fx~, hence 
between 2 and 3 orders of magnitude lower than the total Br. The situation 
is therefore more favourable than in the D — > case. Branching ratios well 
above 10"^ for D° {p° , u") /i'^ ii~ or in the 10"^ range for p^ p~ 

would be indicative of new physics. It is satisfactory to note that present 
experimental bounds are not far above. 

Lastly, we mention the — > 7r"'"'°£+£- decays, whose short distance 
contribution is again related to c — > ui'^£~. In this case, the LD contribution 
reaches [49] a branching ratio of the order of 10"^ in the (^-resonance region 
and a few times 10""^ in the nonresonant region, a situation similar to what 
was encountered in D — > V£^£~ decays. 

5. Be — > -6*7 — a unique opportunity 

The situation described in the previous sections indicates that the prob- 
abihty of observing new physics in D ^ V^, D V£~^£~ or D ^ P£~^£~ is 
rather modest. It would require a mechanism which increases the SD ampli- 
tude of c — > or c — > u£'^£~ by at least one or two orders of magnitude, a 
rather unlikely though not impossible proposition. 
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Fajfer, Prelovsek and Singer [14] have turned to the domain of very rare 
decays and have proposed the idea of exploring the c — > 1*7 transition when 
c is embedded in a beauty particle. In other words, they consider a "beauty- 
conserving" and "charm-changing" decay, which is driven by the c ^ u-y 
transition. As it has been shown by these authors explicitly, such a transition 
has about equal SD and LD contributions, making it an ideal testing ground 
for deviations from SM [14,50]. 

The i?c-iiieson, a compact bound state of two heavy quarks of different 
flavour, c and 6, has been discovered recently at Fermilab [51] and its lifetime 
has been determined as t(5c) = 0.46lo!i6 iO-OSps. The transition c — M + 7 
would lead to the decay — > -|- 7, in which the 6-quark is merely a 
spectator. In order to estimate the SD and the LD contributions to the decay 
one uses the effective Lagrangians of (2) and (4). In (2), the appropriate scale 
for Ct{^) is indeed = rric also for the decay Be i^* -|- 7, and not m^, 
in view of the spectator role of the 6-quark. The general form of the decay 
amplitude is as given in Eq. (5) and we turn now to the calculation of Apc 
and Apy, which have both SD and LD contributions. 
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6. A model for Be -6*7 



The SD contribution calculated from (2) can be expressed in terms of two 
form factors Fi(0), ^2(0): 

e;{B:{p', e') \ma'^-^q^c\Be{p)),.=o = z6'^-^£;4*p>^Fi(0) , (6) 



e;{B:(p', e')\uia''''q,^,c\B,(p)),.=o = [(M|^ - M|.)£* • e * 

-2{e'*-q){p-e*)]F2{0). (7) 

The LD contributions may be separated into two classes related to the two 
terms of (5). The class (I) is related to the a2 term and represents processes 
c — > UQiQi followed by QiQi — > 7, with b as spectator. The qiQi — > 7 transitions 
are expressed by (jiQi hadronization into vector meson, thus we have a vector 
meson dominance (VMD) approximation. The class II of diagrams is related 
to the ai term and corresponds to the quark process cb — > ub with the 
photon attached to quark lines. Only the lowest (pole) states are included 
in the calculation [14,50]. 

The VMD amphtudes of class I are proportional to £* (5* |'U7''( 1—75)0! i?c) 
taken at = 0. This involves one vector and four axial- vector form factors. 
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However, requirements of finiteness at = [46] and gauge invariance imply 
[14] the vanishing of two axial form factors and a relation between the other 
two and accordingly the VMD contribution is expressible in terms of two 
form factors only, V{0) and Ai{0). The amphtudes thus obtained in [14] are 



i.pv — 



+2a2(mc)CvMD 



27r2 
AM 



rric - mu)F2{0) 



Mb, - Met 



(8) 



47r2 



(me + m„)Fi(0) 



+2a2(mc)CvMD 



T/(0) 



Mb, + Mb* 



+ ychV*i,ai{mb) 



X 



I^B.gBiQBz ^ I^BuMbJbJb^ 



Ml, - Ml 



Ml 



Ml 



(9) 



In these expressions the first term is from SD, the second is the LD VMD 
contribution and the third term is the LD pole contribution. Also, 

g'M £(0) 4(0) 



VMD 



2M2 6M2 3M2 

P U) ip 



(-1.2 ± 1.2) X lO-'^GeV 



(10) 



and (y{q.i e)\Vy\Q) = gy{q^)e*^. jii, fi and Qi are couplings related to the axial 
and vector currents and are defined in [14]. 

In order to determine the form factors Ai{0), V{0), -Fi(O), ^2(0) and 
the various iJ>i,fi,gi the authors of Ref. [14] have chosen the nonrelativis- 
tic constituent Isgur-Score-Grinstein-Wise (ISGW) model [52]. This model 



17 



is considered to be reliable for a state composed of two heavy quarks; in 
addition, the velocity of in the rest frame of Be is to a good measure 
nonrelativistic. In the ISGW model the quarks of mass M move under the 
influence of the effective potential V{r) = —Aas/iSr) +c + br with c = —0.81 
GeV, 6 = 0.18 GeV^ [53]. The authors [14] use variational solutions of the 

3 3 —p^r^ 

Schrodinger equation, ip^f) = 7r~ 4/^26 2 for S state with /3 as variational 
parameter. Using accepted values for current quark masses, CKM matrix 
elements and constituent quark masses, one calculates the SD and the LD 
contributions separately, as well as the total branching ratio of Be — > -6*7. 
It is found [14]: 











lol) 




1.7 X K)-'* 








X IC)--' 



Table 3 



As evidenced by the results of Table 3, the SD and LD contributions are 
comparable, which in principle allows one to probe the c — > ^7 transition in 
Be -6*7 decay. Experimental detection of Be -6*7 at a branching ratio 
well above 10^^ would clearly indicate a signal for new physics. It is worth 
mentioning here that at LHC one expects [50] to produce well above lO^Bc 
mesons. 
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Finally, we mention a recent calculation of Aliev and Savci [54] which 
confirms our conclusions [14,50]. They calculate the SD contribution to — > 
B*j by the use of QCD sum rules and find a value for Fj(0) which leads to 
an SD branching ratio for Be — > -6*7 of ~ 1.6 x 10~^, slightly higher than 
presented above, but with the same general conclusions. 

Summary 

We have reviewed the possibility of using various processes to detect de- 
viations from the standard model in the charm sector, using the c — > wy, 
c — > u£'^i~ transitions. The D — > decays are shown to be dominated 
by long distance contributions which usually prevents one from observing 
deviations from the standard model short distance ones. The situation is 
somewhat better in D Vi'^i~ decays, where the gap between SD and 
LD is smaller. Here, branching ratios well above 10"^ for D° — > p°jjL^jjL~ or 
D° — > uj"^'^^" or in the 10~^ range for —>■ p^^"^ ^~ would indicate new 
physics. Of particular interest is the novel decay -8*7 suggested in 

Ref. [14]. In this decay both the SD and LD contributions to the branching 
ratio are in the 10^^ range. The SD contribution is at its natural value. The 
LD one is strongly suppressed, as follows: the Class I VMD contribution 
is very small as a result of the smallness of Cvmd (Eq- 10), which repre- 
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sents a cancellation of vector mesons contributions at a level better than 
10% as a result of GIM and SU(3)i7 symmetry; on the other hand, the class 
II pole contributions is also strongly supressed in view of the appearance of 
the factor VcbV*^, in the cb — > ub pole diagrams. (In D decays we had the 
much bigger KsKTs fa-ctor, which made the LD pole contributions dominant). 
This fortuituous occurrence of SD, LD contributions equality establishes the 
Be — > -6*7 decay mode as an ideal testing ground for physics beyond the 
standard model. To conclude, we stress that this decay has a clear signature: 
the detection requires the observation of a By, decay in coincidence with two 
photons - a high energy one (985 MeV) and a low energy photon (45 MeV) 
in the respective centers of mass of Be and S*. 
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